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Abstract 

We derive formulas for the matrix elements of the two dimensional square lattice Green function along 
the diagonal, and along the coordinate axes. We also give an asymptotic formula for the diagonal elements. 
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I. INTRODUCTION 

In some recent papers HI |2|] Cserti showed how a lattice Green function (LGF) can be used to 
find the resistance between two points in an infinite lattice of resistors. Cserti gives an expression 
for the matrix elements of the LGF in the form of an integral. In this paper we will show how 
to solve this integral for the case of a two dimensional square lattice along the diagonal and the 
coordinate axes. This allows any arbitrary diagonal or coordinate axis LGF matrix element to be 
calculated directly. Formulas for these elements were first derived by McCrea and Whipple 
using a different procedure than that presented here. We will also give an asymptotic formula for 
the diagonal matrix elements that converges to Cserti's asymptotic limit formula for large values 
of n. 

The LGF with which we are concerned here can in general be used to solve the discrete two 
dimensional Poisson equation with boundary conditions at infinity. Therefore it will be useful in 
solving two dimensional electrostatics problems [5] as well as many other problems that can be 
modeled by a Poisson equation. 

n. DIAGONAL MATRIX ELEMENTS 

The matrix elements of the two dimensional square lattice Green function can be expressed in 
terms of an integral as 

f ^ 1 /■'^ /■'^ l-cosn0cosm(^ 

S{n,m) = —^ dQ- (1) 

271^70 Jo 2 — COS 9 — COS (p 



This is essentially the same as Cserti's lU] eq. Bl. We will begin by looking at the diagonal 
elements where m = n. First note that for m = neq. [l]can be rewritten in the following form 

1 /■'^ ,^l-^cos/7((|)-0)-icosn((|) + 0) 



g(n,n) = -^ d(^ dd ^ V — , (2) 

' ^ iTi^Jo ^Jo 2-2cosi((^-0)cosi((^ + 0) 

By symmetry the two cosine terms in the numerator of the integrand can be combined to give 

r ^ 1 r^^r^a l-cosn((^ + 0) 

g(n,n) = —^j dG) dB ^ f (3) 

' ' 2%^ Jo ^Jo 2-2cosi(^-0)cosi((^ + 0) 

In terms of new variables ^' = j{^ — Q),Q' = j{^ + Q) eq. |3lbecomes 

/ ^ l-cos2n0 

g{n,n) = —^ d(^ dB- (4) 

^ ^ An^ Jo ^ Jo l-cos(^cos0 



The integration over (|) can now be done to give 

SV^^^) = irz — — (5) 

An Jo sinG 

Using the identity 1 — cos 2n0 = 2 sin^ «0, eq. Scan be written as 

1 /■'^ sin^nG 1 r,n/^sinn0y 
S{n.n) = — dQ = — / J0 __ sin0 (6) 

271 Jo sin0 271 Jo V sin0 y 

Now if we let x — cos then this becomes the integral of a type II Chebyshev polynomial lO] 

1 /■! 2 1 /"^ 2 

8{n,n) = — U„_i{x)dx=- U„_i{x)dx (7) 
271 J- 1 71 Jo 

The square of a type 11 Chebyshev polynomial can be expressed as 

f/'w = E f/2^w (8) 

k=0 

To prove this identity it is sufficient to show that U^{x) — U'^_^{x) = U2n{x). With x = cos0 we 
have Un{x) = s'm{n + 1)0/ sin and 

sin(n+ l)0-sinn0 cos(n + ^)0 



Unix)-Un-l{x) 



sin0 cos^0 



, , sin(n+ l)0 + sinn0 sin(n + i)0 

Unix) + Un-l{x) = A- = \ , 

sm0 sm50 

rj2(. jj2 _ cos(n + i)0 sin(n + i)0 _ sin(2^ + l)0 _ 

cos^0 sm^0 sm0 
So using eq. [H eq. |7]can be written as 

^(n,n) = i£ Cu2k{x)dx (9) 

Letting x = cos 0, the integrals in eq. |9lbecome 

/ U2k{x)dx= / sin(2A:+l)0t^0 = (10) 

Jo Jo 2A: + 1 

Substituting this into eq. Inland changing the summation index gives 

k=l 
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Note that g{n^n) also obeys a difference equation and that the solution for g{n,n) given in eq. [TT] 
could also be arrived at by solving the difference equation. The difference equation for g{n,n) is 
the same as that given by Cserti |l| eq. 32 for the resistances along the diagonal 

{2n+l)g{n+l,n+l)-Ang{n,n) + {2n-\)g{n~\,n-\) = (12) 

Since the coefficients of this equation add up to zero, if we substitute g{n, n) = Y!k=i /(k), g{n — 
l,n — 1) = g{n,n) — f{n), g{n + l,n + 1) = g{n,n) + 1) into the equation, we will get a first 
order equation for f{n). 

(2n + 1) /(n + 1) - {In - l)f{n) = (13) 

With the initial condition /(I) = ^ the solution to this equation is f{k) = ^ which once again 
gives us eq. [TT]as the solution fox g{n,n). 

We will now derive an asymptotic formula for g{n,n). First note that the g{n,n) elements are 
proportional to the partial sums of a generalized harmonic series. They can also be expressed in 
terms of the standard harmonic series as follows 

.(".") 4 (I I - ill) 4 ("-"H ^''^ 

where we have introduced the notation 

Hn=j^\ (15) 

k=l 

for the nth partial sum of the standard harmonic series. The asymptotic formula for the nth partial 
sum of the harmonic series is [iZJ P- 338 

ft, = ,„„+,+ i__£_||j (16) 

where y = 0.5772156649 ... is the Euler-Mascheroni constant and is a Bernoulli number. Us- 
ing this in eq. [l4lresults in the following asymptotic formula for g{n^ n) 



(n,n) = — 



lnW+y+21n(2) + £ ^^^^^.^ 



52^(22^-1-1: 



(17) 



Without the Bernoulli sum, this is essentially the same as Cserti's [HI] eq. 33 for the asymptotic 
limit of the resistance. 
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III. THE ON-AXIS ELEMENTS 



We now turn to the on-axis elements where m = and eq. [T]becomes 



1 



1 — cosnG 



2 — cos — cos (|) 



(18) 



The integral with respect to (|) can be carried out to give 



1 



dQ 



1 — cosnG 



2nJo v^(2-cose)2-l 



(19) 



Now note that 1 — cosnG = 2sin2(^) and2 — cos0= l+2sin2(|) so that the denominator of the 
integrand in eq. [T9lbecomes ■^J (1 +2sin2(|))2 — 1 — 2sin(|) 1 +sm^{j). Eq. [T^then becomes 



f(n,0) = ^ I dQ 



sin2(f ) 



sin(§)./l+sin2(§) 



Making the change in variable 0' = | , we write eq. I20las 



g{n,0) = l- l^dQ 

Tl Jo 



sinnG 
sinG 



sinG 



l+sin^G 



With X = cos G we then have an integral involving the Chebyshev polynomial Un- 1 (x) 



1 f'Ul,{x) 



dx 



Now if we let x = Vl — cosG then 



1 



g('^,0) = ^y^ t/2_i(Vl-cosG)JG 



1 



2 



^K0) = — £/ t/2KVl-cosG)JG 



If the type I and type n Chebyshev polynomials are expressed in the following forms 

Tn{x) 



1 

2 L 



x+\/x^-\] +[x-\/x^-l 



Unix) 



[x + Vx^y^^ -(^x-s/x^ 



n+l 



then it is straightforward to prove the identity 



2V^23T 



-l)'7^2^+i(v^) 



(20) 



(21) 



(22) 



(23) 
(24) 

(25) 



(26) 
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Using this identity eq. l24lbecomes 

«(".0)4E(-l)V'^^^ii«<*e (27) 

Jo VCOS0 



k=0 

T2k+i (x) can be expressed in series form as 



so that we have 



and eq. |271can be written as 



where 



.5 I itO'-i)" 7 = 2 4 6 

/=1 3 5 7 
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(28) 



= E«(^'i)K;) (30) 



(31) 



Eq. |30|can be used to directly calculate g{n,0) for arbitrary values of n. 
IV. CONCLUSION 

We have derived equations by which g{n,n) and g{n,0) can be calculated for arbitrary values 
of n. For the case of g{n,n) we have an asymptotic formula eq. [iTlthat allows for a quick and 
efficient calculation. In the case of g{n,0) we have eq. l30l whose evaluation can be optimized for 
large values of n. A complete asymptotic formula for g{n,0) has not yet been found. A formula 
very similar to eq. [T7]for g{n,n) has been found for g{n,0) but only the first few terms in the 
Bernoulli sum have so far been determined. Formulas for the general matrix elements g{n,m) 
have been found by us. These formulas are found by solving the partial difference equation for 
g{n,n). This equation can only be solved after a formula for the diagonal elements, g{n,n) has 
been found. 
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